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Decay Rates and Probability Estimates for Massive Dirac 
Particles in the Kerr-Newman Black Hole Geometry 

F. Finster, N. Kamran* J. Smoller| and S.-T. Yau-f 
July 2001 / January 2002 

Abstract 

The Cauchy problem is considered for the massive Dirac equation in the non- 
extreme Kerr-Newman geometry, for smooth initial data with compact support outside 
the event horizon and bounded angular momentum. We prove that the Dirac wave 
function decays in L^^ at least at the rate t~^^^. For generic initial data, this rate 
of decay is sharp. We derive a formula for the probability p that the Dirac particle 
escapes to infinity. For various conditions on the initial data, we show that p = 0, 1 
or < p < 1. The proofs are based on a refined analysis of the Dirac propagator 
constructed in Q. 
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1 Introduction 

The Cauchy problem for the massive Dirac equation in the non-extreme Kerr-Newman 
black hole geometry outside the event horizon was recently studied [^, and it was proved 
that for initial data in L^^ with L^ decay at infinity, the probability for the Dirac particle 
to be located in any compact region of space tends to zero as t — > oo. This result shows 
that the Dirac particle must eventually either disappear into the event horizon or escape 
to infinity. The questions of the likelihood of each of these possibilities and the rates 
of decay of the Dirac wave function in a compact region of space were left open. In the 
present paper, we shall study these questions by means of a detailed analysis of the integral 
representation of the Dirac propagator constructed in Q . This analysis will also give us 
some insight into the physical mechanism which leads to the decay. 

Recall that in Boyer-Lindquist coordinates {t,r,-d,ip) with r>0, 0<'!9<7r, < (p < 
27r, the Kerr-Newman metric is given by Q 

ds'^ = Qjk dx^x 

= ^{dt - a sin2 ^ dipf -u(^ + d^A - ^ (a dt - (r^ + a^) d^f (1.1) 

with 

[/(r,T?) = r'^ + a^cos^^, A(r) = r"^ - 2Mr + a'^ + Q"^ , 

and the electromagnetic potential is 

Aj dx^ = — — - {dt — a sin -d dip) , 

where M, aM and Q denote the mass, the angular momentum and the charge of the black 
hole, respectively. Here a and/or Q are allowed to be zero, so that our results apply also 
to the Kerr, Reissner-Nordstrom, and Schwarzschild solutions. We shall restrict attention 
to the non-extreme case M^ > a? + Q"^, in which case the function A has two distinct 
zeros. 



^•o 



M - Vm2 - a2 - Q2 and n = M + ^J W^ - a^ - Q'^ , 



corresponding to the Cauchy horizon and the event horizon, respectively. We will here 
consider only the region r > r\ outside of the event horizon, and thus A > 0. 

Our starting point is the representation of the Dirac propagator for a Dirac particle of 
mass m and charge e established in y, Thm. 3.6] 

^(t,x) = - Y. due-'"'' ^ t^r ^^"(2^) <^^" I ^o> • (1-2) 

^ fc,nGZ'^~°° a,b=l 

Here ^o is the initial data, uj is the energy, and <.|.> is a positive scalar product (see Q 
for details). The quantum number k arises from the usual separation ~ exp(— z(A; + g)) 
of the angular dependence around the axis of symmetry, whereas n labels the eigenval- 
ues of generalized total angular momentum in Chandrasekhar's separation of the Dirac 
equation into ODEs |Q] . The $^'^"- are solutions of the Dirac equation and t^^" are com- 
plex coefficients; they can all be expressed in terms of the fundamental solutions to these 
ODEs. We postpone the detailed formulas for \I/^'^" and t^^" to later sections, and here 



merely describe those qualitative properties of the wave functions vj/^'^" which are needed 
for understanding our results. Near the event horizon, the ^^'^^ go over asymptotically to 
spherical waves. In the region |a;| > m, the solutions for a = 1 are the incoming waves, i.e. 
asymptotically near the event horizon they are waves moving towards the black hole. Con- 
versely, the solutions for a = 2 are the outgoing waves, which near the event horizon move 
outwards, away from the black hole. Asymptotically near infinity, the ^q*^", |ci;| > m, go 
over to spherical waves. In the region \uj\ < m, however, the fundamental solutions for 
a = 1,2 near the event horizon are both linear combinations of incoming and outgoing 
waves, taken in such a way that ^^'^" and ^l'^" at infinity have exponential decay and 
growth, respectively. 

For technical simplicity, we make the assumption that ^o is smooth and compactly 
supported outside the event horizon. We point out that, while the assumption of compact 
support is physically reasonable at infinity, it is indeed restrictive with respect to the be- 
havior near the event horizon. Furthermore, we shall assume that the angular momentum, 
is bounded in the strict sense that there exist constants ko and no such that 



^(t,x) 



|fc|<fco |n|<no 



dcj e 



-iuit 



E 



t^r *a'^"(x) <*^'^ I *0> 



''ab 



(1.3) 



We expect that the rate of decay is the same if an infinite number of angular modes are 
present. Namely, away from the event horizon, modes with large angular momentum feel 
strong centrifugal forces and should therefore be quickly driven out to infinity, whereas 
the behavior near the event horizon is independent of the angular momentum. However, 
it seems a very delicate problem to rigorously establish decay rates without the assump- 
tion (|1.3| ) , because this would make it necessary to control the dependence of our estimates 
on k and n. Finally, we assume that the charge of the black hole is so small that the grav- 
itational attraction is the dominant force at a large distance from the black hole. More 
precisely, we shall assume throughout this paper that 

mM > \eQ\ . (1.4) 

We now state our main results and discuss them afterwards. 

Theorem 1.1 (Decay Rates) Consider the Cauchy problem 

{i-iWj-m)^{t,x) = 0, ^(0,a;) = ^o{x) 

for the Dirac equation in the non-extreme Kerr-Newman black hole geometry with small 
charge ^l-^j )- Assume that the Cauchy data ^q "is smooth with compact support outside 
the event horizon and has bounded angular momentum (l-S). 



(i) If for any k and n, 

limsup|<^^'^"|^o>| / 

uj'\m 

then for large t, 

\^{t,x)\ 

with c = c{x) 7^ and any e < ^. 



or 



limsup|<^^'^"|^o>| / 0, (1.5) 

u}/'—m 



ct-6 +0{t- 



(1.6) 



(ii) If for all k, n and a = 1,2, 

'0 



<^^'^" I ^n> = 



for all u in a neighborhood of zizm, then |^(t, x)| has rapid decay in t (i.e. for any 
fixed X, ^{t,x) decays in t faster than polynomially) . 

Theorem 1.2 (Probability Estimates) Consider the Cauchy problem as in Theorem 



1.1, with initial data ^o normalized by <^q\'^q> = 1. Let p be the probability for the 



Dirac particle to escape to infinity, defined for any R > ri by 

p = lim (^j^^){t,x)uj dfi , (1.7) 

*-"°° J{R<r<oo} 

where v denotes the future directed normal to the hypersurface t = const and dfi denotes 
the induced invariant measure on that hypersurface. Then p is given by 

p = - Y] E / d^f^-2|t^ri'l <*2""i^o>' . (1.8) 

Accordingly, 1 — p gives the probability that the Dirac particles disappears into the event 
horizon. Furthermore, 

(i) Suppose that the outgoing energy distribution for \lo\ > m is 

non-zero, i.e. 

<^fc'^" I ^o> / 

for some lo with \uj\ > m. Then p > 0. 

(ii) If the energy distribution of the Cauchy data has a non-zero contribution in the in- 
terval [—m,m], then p < 1. 

(iii) // the energy distribution of the Cauchy data is supported in [—m,m], then p = 0. 



(iv) // (1.5) holds, then < p <1. 



The decay rate of t~6 obtained in Theorem |1.1| quantifies the effect of the black hole's 
gravitational attraction on the long-time behavior of massive Dirac particles. Before dis- 
cussing this effect in detail, it is instructive to recall the derivation of the decay rates 
in Minkowski space. We denote the plane-wave solutions of the Dirac equation by ^r , 
where k is momentum, e = ±1 is the sign of energy, and s = it refers to the two spin 
orientations. The plane-wave solutions are normalized according to 

(^fc.el%sv) = S{k-k')6,,,6,,,, 
where (.|.) is the usual spatial scalar product 

The Dirac propagator is obtained by decomposing the initial data into the plane-wave 
solutions, 

*(t,f) = Y. I dk ^kse(t,S) (^fcs.(^ = 0) I ^O) , 



and a straightforward calculation using the explicit form of the plane- wave solutions yields 
that 

^{t, f ) = 2vr I ^ (^ + m) 6{e - m") r(a;) e^^^^ 7° ^o(^) , (1.9) 

where ^o(^) is the Fourier transform of '^q{x) (and as usual x = (t,x), k = {uj,k), 
1^ = kj^^ , and T is the step function T{x) = sgn(x)). Let us assume for simplicity that 
the initial data is a Schwartz function. We write ( |1.9D as a Fourier integral in lo, 

/oo 
^iuj,x)e'''^^duj (1.10) 

-00 

where 

—-^ (^ + m) 5{u;^ - \k\^ - m") r(a;) e^^^ 7° ^o(fc) • 

We consider the w-dependence of ^ for fixed x. The (5-distribution gives a contribution to 
the momentum integral only for k on the sphere \k\'^ = uP' — rr? . Thus ^(o;, x) vanishes for 
|u;| < m and has rapid decay at infinity. Furthermore, ^ is clearly smooth in the region 
|u;| > m. For \uj\ near m, ^ has the expansion 



4^2 



87r2 



(w + m7°) ^o(O) \/u;2 - m2 + ©(^^ - m^) . 



A typical plot of |^(u;,x)| is shown in Figure |^(a). If ^0 vanishes in a neighborhood of 
A; = 0, then ^(.,x) is a Schwartz function, and thus its Fourier transform ( [I.IOD has rapid 



decay. This is the analogue of Case (ii) of Theorem LI. However, if ^o(O) 7^ 0, the decay 
rate is determined by the square root behavior of ^ for |lij| near ra. A change of variables 
gives that for any test function t] which is supported in a neighborhood of the origin. 

An integration-by-parts argument shows that the last integral is bounded uniformly in t, 
and is non-zero for large t if 77(0) 7^ 0. From this we conclude that in Minkowski space, 

3 

|^(i, x)| decays polynomially at the rate i~2. 

We now proceed with a more detailed discussion of our results, beginning with the rates 
of decay obtained in Theorem |l.l| . Naively speaking, a massive Dirac particle behaves near 
the event horizon similar to a massless particle, i.e. like a solution of the wave equation. 
In Minkowski space, solutions of the wave equation decay rapidly in time according to the 
Huygens principle. On the other hand, at large distance from the black hole the solutions 
should behave like those of the massive Dirac equation in Minkowski space, which decay at 
the rate i~ 2 . It is thus tempting to expect that the solutions of the massive Dirac equation 
in the Kerr-Newman black hole geometry should decay at a rate which "interpolates" 
between the behavior of a massive particle in Minkowski space and that of a massless 



particle, and should thus decay at a rate no slower than t~2. However, Theorem 1.1 
shows that this naive picture is incorrect, since the rate of decay we have established 
for a massive Dirac particle in the Kerr-Newman black hole geometry is actually slower 
than that of a massive particle in Minkowski space. Thus the gravitational field of the 



1951 





Figure 1: Typical plot for ^ in Minkowski space (a) and in the Kerr-Newman black hole 
geometry (b). 



black hole affects the behavior of massive Dirac particles in a more subtle way. One can 
understand this fact by comparing the plots in Figure |l], which give typical examples for 
the energy distribution of the Dirac wave function in Minkowski space and in the Kerr- 
Newman geometry. One sees that in the Kerr-Newman geometry, there is a contribution 
to the energy distribution for |ti;| < m, which oscillates infinitely fast as cv approaches m. 

_ 5 

When taking the Fourier transform, these oscillations lead to the decay rate t e given 
Theorem ^J (see the rigorous saddle point argument in Lemma |3.3| ). 

The oscillations in the energy distribution in Figure 0(b) are a consequence of the field 
behavior near spatial infinity. On a qualitative level, they can already be understood in 
Newtonian gravity and the semi-classical approximation. Namely, in the Newtonian limit 
of General Relativity, the momentum A; of a relativistic particle is related to its energy uj 

by 



H 



UJ + 



■) 



m 



Thus the particle has positive momentum even if uj < m, provided that the Newtonian 
potential is large enough, ^^^ > m — uj. This means in the semi-classical approximation 
that the wave function ^(r) has an oscillatory behavior near the black hole, 



^(r) ~ exp I ibz 



k ds 



for 



r < R = 



niM 



m — ijj 



and will fall off exponentially for r > R. As a consequence, the fundamental solutions 
•(^kum £qj, 1^1 ^ j^ involve phase factors ~ exp(ibi J kds). In the limit co y m, R ^ oo, 
leading to infinitely fast oscillations in our integral representation. This simple argument 
even gives the correct quantitative behavior of the phases ~ (m — io)~2 . 

The fact that the decay rate in the presence of a black hole is slower than in Minkowski 
space has the following direct physical interpretation. One can view the gravitational 
attraction of the black hole and the tendency of quantum mechanical wave functions to 
spread out in space as competing with each other over time. The component of the wave 
function for uj near m and uj < m, which is responsible for the decay rate t~e, has not 
enough energy to propagate out to infinity. But since it is an outgoing wave near the 
event horizon (note that in (|1.5|) the fundamental solutions \I/^'^" enter only for a = 2), 
it is driven outwards and resists the gravitational attraction for a long time before it will 
eventually be drawn into the black hole. As a result, the Dirac particle stays in any 



compact region of space longer than it would in Minkowski space, and thus the rate of 
decay of the wave function is slower. 

According to this interpretation, our decay rates are a consequence of the far-field 
behavior of the black hole. Similar to the "power law tails" in the massless case (see |p), 
our effect can be understood as a "backscattering" of the outgoing wave from the long- 
range potential, but clearly the rest mass drastically changes the behavior of the wave 
near infinity. We expect that result for the decay rates should be valid even in a more 
general setting, independent of the details of the local geometry near the event horizon. 
Furthermore, the decay rates should be independent of the spin. This view is supported 
by P, 0], who obtained the rate t~6 for massive scalar fields in a spherically symmetric 
geometry using asymptotic expansions of the Green's functions. 



Theorem |1.2| gives a precise formula for the probability that the Dirac particle either 
disappears into the black hole or escapes to infinity. In cases (i)-(iv) we give sufficient 
conditions for these probabilities to occur. These results are consistent with the general 
behavior of quantum mechanical particles in the presence of a potential barrier and can be 
thought of as a tunnelling effect. In case (iii), the particle does not have enough energy to 
escape to infinity. Thinking again in terms of a tunnelling effect, the Dirac particle cannot 
tunnel to infinity because the potential barrier (which has finite height m) has infinite 
width. Finally, one might ask whether p = 1 can occur; i.e., that the particles escapes to 
infinity with probability one. This is indeed the case for very special initial data, whose 
energy distribution is supported outside the interval [—m,m\ (see Corollary |9.3| below). 

We conclude by remarking that a number of significant results are known for the long- 
term behavior of massless fields in black hole geometries. These results do not capture our 
effect, which is intimately related to the presence of a mass gap in the energy spectrum. 
Price 1^] discussed the rates of decay of massless fields in the Schwarzschild background 
for special choices of initial data. His decay rates depend on the angular momentum and 
are faster than the ones we have derived. A rigorous proof of the boundedness of the 
solutions of the wave equation in the Schwarzschild geometry has been given by Kay and 
Wald [^. Beyer pursues an approach using C*^-semigroup theory, which also applies to 
the Kerr metric and the massive case |l|. An important contribution to the long-time 
behavior of gravitational perturbations of the Kerr metric has been given by Whiting Q. 

2 The Long-Time Dynamics under a Spectral Condition 

We begin the analysis with the case when the energy distribution of the Cauchy data is 
zero in a neighborhood of a; = ztm. The following theorem is an equivalent formulation 
of Theorem |l.l| (ii). 

Theorem 2.1 Consider the Cauchy problem 

{i-fWj-m)'^{t,x) = , ^(0,x) = ^0(2;) 

for smooth initial data with compact support outside the event horizon. Assume that 
angular momentum is hounded and that the energy is supported away from uo = ibm, i.e. 

^0 = - E E / +/ +/ M^ E*ar^a""<^^"I^O> (2.1) 

|fe|<fco \n\<nQ a, 0=1 

for suitable e > 0. Then for all x, ^{t,x) has rapid decay in t. 



Before giving the proof, we recall a few basic formulas from [Q] . The separation ansatz 

kum 



for the fundamental solutions '^'^"- is 



*^-"(t,r,^,V.) 



-iuit p-i(k+^)ip 






(2.2) 



where X = (X+,X_) and Y = (y+,y_) are the radial and angular components, respec- 
tively. The radial part X{u) is a solution of the radial Dirac equation M, eqn (3.7)] 



d^+^^(^Ho -1 



X 



J.2 _j_ g2 



imr — A 

-imr — A 



X 



(2.3) 



where 



0(u) 



UJ + 



{k + \)a + eQr 



J.2 _j_ Q,2 



A 



2Mr + a^ + Q^ 



A is the angular eigenvalue (which depends smoothly on w), and u G (— cxd,oo) is related 
to the radius by 



du r^ + a^ 



dr A 

To analyze X in the asymptotic region u ^ —oo, one employs for X the ansatz 



X{u) 
and obtains for / the equation 

i{QQ-n{u)) 



^i^QU 



r{u) 



(2.4) 



(2.5) 



du 



1 
-1 



+ 



A 



J.2 _|_ q2 







-2iQu 



_2if2n 



{—imr — A) 



(imr — A) 




/ 



(2.6) 



Standard Gronwall estimates yield that the fundamental solutions of ( |2.3| ) have the asymp- 
totic form y, Lemma 3.1] 



XJu) 






+ Ro{u) , 



where |i?o(^)| < c exp(du) for suitable constants c, d > and 



/( 



01 



/o; 



(2.7) 



(2.i 



In the asymptotic region u ^ oo, one transforms the spinor basis with a matrix -B(u) such 
that the matrix potential in ( |2.3D becomes equal to the diagonal matrix —iQ,{u)a'^ {a^ are 
the Pauli matrices). One must distinguish between the two cases |u;| < m and |a;| > m. In 
the first case, i^{u) is imaginary for large u, and thus there are two fundamental solutions 
Xi and X2 with exponential decay and growth, respectively, and we normalize them such 
that 

lim \X{u)\ = 1 . (2.9) 

u— >— 00 



In the case \uj\ > m, 0,{u) is real for all u. The ansatz 

X = B (^ ''J^f^^^ ) with <J>'(n) = n{u) (2.10) 



gives the differential equation 



— / = M(u) f with \M(u)\ < - 

du u 



—f = M{u) f with |M(n)| < ^ , (2.11) 



which can again be controlled by Gronwall estimates. Thus one obtains the asymptotic 
formula Q, Lemma 3.5] 

/ cosh 6 sinh0 \ / e~**(") f+ \ 

where |-Roo| < c/u for suitable c > and 

e = - log , ^ = T{uj)(^uj^-m^u+ 7 logu . 2.13 

The complex factors /^ ^ in ( |2.12| ) are the so-called transmission coefficients. Furthermore, 
we introduced the functions ta{a), < a < 2tt, in terms of the transmission coefficients 
by i eqn (3.47)] 

h{a) = /+2e'*" - f^2e"' , h{a) = -/+ie— + /" ^ e^ . (2.14) 

Finally, the coefficients (tafe)a,6=i,2 are given by 



tab 



6ai 6bi if \lo\ < m 

1 f- t,h , .,, , (2-15) 

da II \uj\ > m. 



[ 2tt Jo |tiP + |t2P 



Proof of Theorem 2.1. Since ( p.l| ) contains only finite sums, we can fix k,n and consider 
one summand. The coefficients in the differential equation (|2.6| ) are smooth in w, and their 
^-derivatives are integrable on the half-lines u G (— oo, uq] for uq sufficiently small. Apart 



from the singularities at u; = ±771, the same is true for the differential equation (|2.11 ) for 



u on the half line [ui,oo) and ui sufficiently large. Since furthermore the ansatz ( |2.5| ) 
is smooth in u; G M \ ((— m — e, —m + e) U {m — £,m -\- e)) and ( ^.lOD is smooth in w G 
(— 00, — 77Z— e]U[?7i+e,oo), and also the coefficients of the differential equation ( |2.3| ) depend 
smoothly on co in the bounded interval u £ [uo,ui], we conclude that the fundamental 
solutions ^^°'^"°(x) and the transmission coefficients fa"'° depend smoothly on a; G 
M\{{—m — e,—m + e)L){m — e,m + e)). Hence the integrand in ( |2.1D is a smooth function 
in uj (which vanishes for |a;| — m < e). Since ^0 has compact support and the fundamental 
solutions for a; ^ cxd go over to plane waves, it is clear that the w-derivatives of the 
integrand in ( p.l[) are all integrable. It follows that the Fourier integral (|1.2| ) has rapid 
decay. ■ 

According to this theorem and using linearity, it remains to analyze the energy distribution 
in a neighborhood of a; = ±771. Since all constructions and estimates are similar for positive 
and negative co, we can in what follows restrict attention to a neighborhood of u; = +m. 



3 Decay Rates of Fourier Transforms — Basic Considerations 



In this section, we derive estimates of some elementary Fourier integrals. Our decay rate 



t 6 ultimately comes from Lemma 3^. We always denote by e a parameter in the range 
< e < 4. 
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Lemma 3.1 Let g G L°°{M.) n C^((0,oo)) with compact support and assume that for a 

suitable constant C , 

C 

\g'{a)\ < — for alia > 0. (3.1) 



a 



Then there is a constant c = c{g) such that for all t > 0, 



/■oo 

/ e*"* g{a) da 
Jo 



< Ct'6^ 



Proof. Assume that supp g C [—L, L]. For given 5 > 0, we split up the integral as 



e*"' g{a) da 







e*"* g{a) da + / e*"* g{a) da . 
Js 



The first term can be estimated by 



1-5 

/ e'''^g{a)da 
Jo 



< ci6 



with ci = sup \g\. In the second term, we integrate by parts, 

1 f°° / d 



I 



e'"* g(a) da = - — e'"* g(a) da 

S it Js \da ) 

= --e''g{8) - - / e-V(«)d«, 
rt it Js 



and estimate using (|3.1|), 



e*°* g{a) da 



< J + y(logL-log5). 



We choose 6 = t e ^ to conclude that 



e*"* g{a) da 



ci C 
< cit"^"^ + — H (logL-logt^^^^) 



and this has the required decay properties in t. ■ 

In the next lemma we insert into the Fourier integral a phase factor which oscillates 
infinitely fast as a \ 0. 



Lemma 3.2 Let g be as in Lemma 3A. Then there is a constant c = c{g) such that for 
all t > 0, 



exp I iat j= g(a) da 

V \/a 



< Ct~6- 



(3.2) 



10 



Proof. We set 



Then 



</>(«) 



,// N Is 



at 



a 



4>"{a) 



3 _5 

-— a 2 
4 



(3.3) 



We integrate the Fourier integrals by parts, 



s*'^(") g{a) da 






t e 

'o 






) 9__9: 



kl2 



da , 



and obtain the bound 



e^'^(") g{a) da 



< 



g' 9<P" 



A'2 



da 



{\9'\\<P'\--^^' + \9\\r\m-'^^'}da . 



< t 6 ^ / ( ci 15 I a4 2^ + C2 l^l a 4 i^Xda 



According to (^.31), we can estimate the factor \^'\ e ^ from above by t s ^, whereas for 
the factors \(\)'\ in the curly brackets we use the bound \^'\ > ^ a~2. Furthermore, we 
substitute in the formula for (p" in ( |3.3[ ) and obtain 



s*'^^") 5(a) da 



with two constants ci and C2. Using that g is in L°° and g' satisfies the bound ( |3.1| ), one 
sees that the pole in the last integrand is integrable. ■ 

The following lemma deals with the Fourier integral when we replace the minus sign in 
the integrand of (3.2) by a plus sign. Reversing this sign completely changes the long-time 
asymptotics. We estimate the Fourier integral using a rigorous version of the "saddle point 
method." 



Lemma 3.3 Let g be as in Lemma p. j| . Then there are constants c = c{g) and ci 

4 1 

2~3 3^2 y^ such that for all sufficiently large t, 



exp ( iat -\ — "— ) g{a) da — ci e**^" 5(00) t 
V va. 



< Ct"6" 



where ao and 0o are given by 

ao = (2t)-i , 
Proof. We introduce the function (p by 



4>o 



9 \ 3 

4* 



(3.4) 



(3.5) 



a 



at + 



a 



Then 



,,, , 1 _3 

(j) [a) = t — - a 2 



<P"ia) 



3 _5 

— a 2 
4 



(3.6) 



11 



One sees that 4>{a) has a minimum at oq with (j){aQ 

ct)"{ao) 
We set 



and 



1 5 

2-3 3t3 . 



(3.7) 



t"6 



+S 



For large t, 5 < oq. We spht the integration range into two regions Di and D2 with 

Di = (0, ao - 6)U (ao + 5, 00) , -D2 = ["o - S,aQ + 6] . 

Let us first estimate the integral over Di. An integration- by-parts argument similar to 
that in the proof of Lemma pl3 gives 



e^"^ g da 



Di 



< 



g{ao + S) 



(t>'{ao + 6) 



+ 



g{ao - 6) 



<p'{ao - 6) 



+ 



D, 



g' 9<\>" 



m 



da . 



Putting in the formulas for (p' and (p" given in ( |3.6D , and using that g £ L°^ together 
with (^j|), one sees that for suitable c, 



e^"^ g da 



Di 



< Ct~6' 



Next we show that the leading contribution to the integral over D2 is given by the saddle 
point approximation. To this end, we introduce the quadratic polynomial 

(f)s{a) = (f)Q + - (j)"{ao) (a - oq)^ . 
Then the mean value theorem gives for sufficiently large t, 

e^'P g - e'^^ g{ao)) da < snp{\ g' \ + \{cl> - <Psy g\) 6^ 



ao 



where in the last step we used that 5e < g. Finally, we compute the contribution of the 
saddle point approximation. 

/ e'^Sg(ao)da = e'^'^ giao) [ e^ "^"("o) "' da . 
Jd2 J-s 

Litroducing the new variable s = 2 0"(oo) ct^ gives 



e*<^« g{ao) da = e*<^« g{ao) 



i4>0 



D2 



r{ao) 



ds 



with 



L = -(l)"{ao)S^ = 2-1 3 1^' 



(3. 



Using ( |3.7| ), we conclude that 

e^'^« g{ao) da = e''^° g{ao) ^1 {a + R{t)) 



D2 
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with 



Jo Vs 
Rit) 

The error term R{t) can be integrated by parts, 

ds 



-23 3~2 / -—ds . 
Il VS 



5 1 



1 1 /■ _3 , 

< ^= + - / S 2 (is _ 



and this shows according to ( |3.8D that R{t) decays in t at the desired rate, \R{t)\ < ct~^. 



4 The Planar Equation 

Let us transform the radial Dirac equation (|2.3| ) into an equation for a real 2-spinor as 
follows. We first unitarily transform the spinor X according to 



X ^ X = UX with 

Then X satisfies the equation 



U 



P 3 
expl t-a-^ 



P 



arctan ■ 



A 



mr 



du 



with 



a(u) 
b{u) 



n{u) + 



-a b 
-b a 



Xm 



X 



A 



m^r^ _l_ ^2 2(r2 + a'^) 



J.2 _|_ Q^ 



y rn^r^ + A^ . 



(4.1) 

(4.2) 

(4.3) 
(4.4) 



Notice that the transformation \J is regular for all n G M, and that the second summand 
in (|4.3|) has nice decay properties for u — > ±oo. Next we employ the ansatz 



X 



Via) ^+ - i^- 
-ip"^ — i r(a) ^" 



i; 



1 / T{a) X+ - X- 
2\iX+ + iT{a)X- 



with a complex 2-spinor %[). Then ip satisfies the equation 



d_ 
du 



^ 



-g 
f 



^ 



with 



/ 



|a| + 6 , 



9 



(4.5) 



(4.6) 



(4.7) 



The coefficients in (|4.6D are all real, and so we can study the real and imaginary parts of 
ijj separately. Thus we assume in what follows that tp is real and then call (^^) the planar 
equation. 
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We bring the planar equation into a form more appropriate for our estimates. For 
given uq we introduce the new variable 



x{u) 



VlMi^dT 



and set 



log 



In the case g > 0, (|4.6| ) becomes 



ax 



-e 




^l^ 



Employing the ansatz 



iP 



e 2 



h 

6 2 COS 

_h 

e 2 sin 



2 

x+i} 
2 



with real functions L{x) and 4>{x) gives the equation 

-L' ai + h'b + {1 + ^') 02 



0.2 



with 



ai 




sm 



X + S 



02 



COS^ 



cos 
— sin 



.-, x+i 



2 

x+S 
2 



Elementary trigonometry shows that 

b = cos(x + •&) ai — sin(x + t?) 02 . 
Hence the planar equation takes the form 



1}' = h' sin(x + {}) 
In the case 5 < 0, the ansatz 



L' = h' cos(x + t?) 



V' 



e 2 



6 2 cosh^i^ 



6 2 sinh 



x+i? 



gives similarly the equations 

{} = h sinh(x + •&) 



L = h cosh(j; + •&) 



(4.8) 



(4.9) 



(4.10) 



(4.11) 



(4.12) 



(4.13) 



(4.14) 



We can now give the strategy for the proof of Theorem |l.l| (i) . First, in the next section, 
we will obtain estimates which will enable us to control the function h' which appears in 
the planar equations ( [4.12| ) and ( 4.14 ). Then we will carefully analyze the solutions {^,L) 
of these planar equations, and this will allow us to study the time-dependence of the 



propagator (1.3). For the analysis of the planar equations, it is necessary to consider both 
cases to > m and to < m separately; this will be done in Sections ^ and |^, respectively. 
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5 Uniform Bounds for the Potentials 

In this section, we shall derive estimates for the function h{x) (as introduced in ( [4.9D with 
X according to ( |4.8D ) as well as for its partial derivatives with respect to x and to. The 
usefulness of our estimates lies in the fact that they are uniform in to for w in a small 
neighborhood of m, oj £ (m — 5,m + 5). The main technical difficulty is that x is defined 
via an integral transformation ( |4.8| ), and thus h{x) depends on a; in a nonlocal way. On the 
other hand, our estimates also show the advantage of working with the variable x. Namely, 
by introducing x, the w-dependence of h becomes small in the critical regions near infinity 
and near the poles of h, in the sense that (a; — m)d^h[x) has bounded total variation in 
a;, uniformly in uo. This will be essential for getting control of the w-dependence of the 
solutions to the planar equation (see Lemmas p^ and |7.5|) . Since the technical details of 



the proofs of Lemmas 5T, 5^, and ^^ will not be needed later on, the reader may consider 
skipping these proofs in a first reading. 

In what follows, we often denote derivatives by a lower index, e.g. h^ = d^h. Further- 
more, we denote constants which are independent of ui by c; the value of c may change 
throughout our calculations. For clarity, we sometimes add a subscript to c to mean a 



fixed constant. According to their definition ( |4.7| ) and ( |4.3| ),(4.4), the functions / and g 
have for large u the expansion 

mM -eQ ,^( \ 

u \n^ 

, , mM + eQ ^^( I 

Our notation 0{u~^) implies that the error terms depend smoothly on w, and that their 
^-derivatives have the natural scaling behavior, i.e. 

5^0(n-") = ©(-u-") and 5„0('u-") = ©(n^^-^) . 

Our assumption ( |1.4D ensures that for large n, g is monotone decreasing, whereas / is 
increasing. 

We begin with the case uj > m. In this parameter range, we fix uq independent of lo. 
By choosing uq sufficiently large, we can arrange that the following estimates hold. 

Lemma 5.1 There are constants c,6>0 such that for all to € (m, m + 5) and x > 0, 

< -h'{x) < ^— (5.1) 

^ ' ~ l+x ^ ' 

|/i"(x)| < 7:^^ (5.2) 



/o UJ — ra 

Proof. We set e = to'^ — m? and introduce the function 

Then h and p have the asymptotic expansions 



(l + x) 
\h'^{x)\ dx < — - — . (5.3) 



h{u) = ^logf^ ^-^ + - + o(\\] (5.4) 




p{u) = 2.e + ^ + 0[^\ (5.5) 
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with positive constants a and /?, which depend smoothly on to and are bounded away from 
zero as e — > 0. Our first step is to bound the function x{u), ( |4.8| ), as well as its inverse 
u{x). According to ( ^.51) , there are (possibly after increasing uq), constants 01,02 > 
such that 

2We + — < p{u) < 2,/e + — (5.6) 

\ u \ u 

for all CO G {m, m + 6) and u > uq. We introduce the functions x and x by 
X = 2\Ju (a\ + eu) — bi , x = 4y^u (02 + eu) — 62 , 

where the constants 61 and 62 are chosen such that x{uq) = = x{uq), 



h = 2v^uo (ai +euo) , 62 = W'^o ("2 + e^^o) • (5-7) 



Then 



x{u) = , . , . < 2W < p{u) 

yju{ai +eu) V u 

x{u) = , , , , > 2W > p{u) , 

and integration yields that x and x are bounds for x, 

x{u) < x{u) < x{u) for all u > uq. (5-8) 

The functions x and x are strictly monotone and thus invertible. Their inverses are 
computed as follows, 

xT^ix) = — (ya? + e(x + 6i)2 -ai 

1 e{x + bif 1 (x + 6i)2 



< 



2e ^ai + e{x + bi)^ + ai ~ 2 ^ {x + bi) + ai 



X (x) = 



^ (^J44 + e (x + 62)2 - 202) 



1 e(x + 62)2 > 1 (x + 62)2 



4e v/4ai+e(x + 62)2 + 2a2 4 ^e (2; + 62) + 4a2 ' 

where in the last step we applied the inequality V« + b < ^/a + Vb {a,b > 0). The 
inequalities (|5.8| ) yield for the inverses that x~^(x) > u{x) > x~^(x). Thus the functions 
u and u defined by 

M = i ^" + ^^)' u(x) = i ^^ + ^^)' (5 9) 

-^^^ 4 V^ (x + 62) + 4a2 ' ^ ' 2^{x + b,) + a, ^^'^^ 

are bounds for tt(x), 

t((x) < u{x) < n(x) for all x > 0. (5.10) 
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Let us derive (^.ip. Since / and g are monotone increasing and decreasing, respectively, 
h{u) is decreasing and thus h'{x) < 0. Furthermore, 



\h'ix)\ 



- \h'{u)\ = ^ 

P 2/9 



g'iu) fiu) 



9 



f 



-^ Wiu) f - f'{u) g\ < -^-^ 



(5.11) 



We employ (5^), ( 5.1C| ), and (5^) to obtain 



\h'{x)\ 



< 



— u 2 (^ai + eu) 



1 Cl _i , , 

2 < — u 2 (ai + eu) 



< 



2ci 
2ci 



^/e{x + 62) + 4a2 
(x + 62)^ 



\/e {x + 62) + 4a2 



e(a; + 62)^ + 4ai ^(x + 62) + 16aia2 



X-\-h2 



{y/e {x + 62) + 4a2)^ 
(e(x + 62)^ + 160102) 



The square bracket is bounded uniformly in e and x, proving (|5.l| ). The second derivative 
of h is computed to be 



h"{x) 



-^(^(9'{u)f-f'{u)g)) 
p du \p'^ J 



-A (p'(n) (5'(n) / - /'(n) 5)) + -^ j^{g'{u) f - fiu) g) . 



Thus 



\h"{x)\ < 



C2 



u'^p^ 



+ 



C2 



«3p4 



(5.12) 



and ( |5.6D shows that 

|/i"(x)| < ^C2U'^ {ai+eu)'" + C2U 
We substitute in (|5.9|), 



-W. , .„,A-3 , ._„,-! („^+^^)-2 < ._„,-! 



C3M 



ai + eu) 



\h"{x)\ < 



I6C3 



(xi + 6i)2 



\3 H 



(V^(x + 62) + 4a2 
(e (x + 62)2 + 16aia2)2 



The square bracket is again bounded uniformly in e and x, and this gives (|5.2|) . We finally 
estimate h'^{x). Since the relation between w and e is one-to-one and smooth, we can just 
as well consider the e-derivative h'^{x). Since h{x) is not given in closed form, we need to 
compute /ie(2;) via the formula 



eh^{x) = eh^{s) + eh'[s) s^{x) , 



(5.13) 



where s = s{u) is a suitable variable. Clearly, hs{x) is independent of how s is chosen. 
However, if we take for s too simple a function (e.g. s = u), then it turns out that 
h^{s) will develop singularities in the limit e ^ 0, which are compensated in ( |5.13| ) by 
corresponding singular contributions to the second summand, making the analysis very 
delicate. To bypass these difficulties, it is convenient to choose for s{u) a function with 
a similar qualitative behavior as x{u); this will make it possible to estimate the two 



summands in (|5.13D separately. We set s = [li-inverse + 61], so that u = u{s{u) — bi); i.e. 

h{s) = h{u = u{s — bi)) 
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with u according to (5^) and s E [61, cci). Then the expansion ( |5.4D becomes 
Ks) = \ log [e Ai + ^^ + ^ + ^/i 0{s-^) + 0(s-3)) 



with positive constants Aj which depend smoothly on ^fe and are uniformly bounded away 
from zero. Differentiating with respect to e gives 

. . ^ 1 Ai£g^ + A2 ^^5 + ^^0(5°) 

'^""^ ~ 2 Aies2 + 2A2^/^s + A3 + ^/^O(s0) + O(s-i) ■ 

We want to show that this function has bounded total variation. To this end, we differ- 
entiate with respect to s and obtain 



eh!,{s) 



ci eis^ + C2 es + C3 ^/e + et 0{s) + e 0{s^) + ^ ©(s"^) 



2 (Ai es2 + 2A2 \/Fs + A3 + ^/F O(s0) + 0{s-^)Y 
Hence by choosing b small enough and uq (and thus 62) large enough, we can arrange that 

I, /I ^ e2s^ + es + ^ 

P \n < c — 

' ^' - (Aie.2 + A3)2 

for all e and all s > 61. The L^ norm of the rhs is bounded uniformly in lo. Namely, 
setting t = VAie s shows that for n = 0, 1, 2, 



iVesY 



b2 



71 + 1 /"oo 



(Aies2 + A3)2V - 1 J^ (t2 + A3)2 



(it 



and the last integral is finite, independent of e. It remains to estimate the total variation 
of the second summand in ( ^.13 ). More precisely, in order to finish the proof of ( |5.3| ), we 
shall show that 

^{h'{s)ese{x)) 



dx 



dx < c . 



(5.14) 



We first derive sufficient conditions for ( ^.14| ). The relations 

= deS{x{s)) = Se{x)+s'{x)Xe{s) h' {x) = h' {s) s' {x) 

yield that 

h'{s)Se{x) = -h'{x)Xe{s) . 

Differentiating with respect to x, one sees that it suffices to bound the L^ norms of the 
expressions 

and h"{x) e 2;e(s) 



h'{x) '-4^ 

x'[s) 



uniformly in e. Substituting in the bounds (5.1) and (^.2|), we conclude that the following 
inequalities imply that (5.14) holds, 



1 // ^ 

- < x'[s) < c 

'^e\x',{s)\ 







l + x 

£ \Xe{s)\ 
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dx < c 



dx < c . 



(5.15) 
(5.16) 
(5.17) 



We begin the proof of ( ^■15|) -( ^.17 ) by computing x'{s 

A short calculation using ( |5.9| ) and ( |5.5| ) gives 

s{^/es + 2ai 



p{v) dv = p{u{s — bi)) if {s — bi) 



(5.18) 



u'{s - 61) 
pHuis-b,)) 



2(VFs + ai)2 



> 



1 



(4es2 + 8/3^s + 8/3ai) + ^ 0{s-^) + 0(5" 



and thus 



x'is? 



(Vis)^ + gg (^/is)^ + • • • + gp + gi ^(g^) + e 0{s) + ^i 0(^O) + 0{s~^] 
( VFs)4 + 63 (^s)3 + • • • + 60 



(5.19) 

with coefficients aj,bj > and go, 60 > 0. Possibly after increasing uq and decreasing 
5, we can neglect the error terms. The fraction in (5.19) is clearly uniformly bounded 
from above and below. This proves ( 5.15| ). Integrating ( 5.15 ), we obtain that the ratio 
(1 + x)/s is uniformly bounded from above and below, and thus in ( 5.16| ) and ( ^.17 ) we 
may replace the factors (l + x) by s. Using ( ^.151 ) we may furthermore replace the integral 
over X S (0, 00) by the integral over s G (61, 00). Next we differentiate ( |5.19| ) with respect 
to e. A short computation shows that 



ex',{s)x'{s) 



_a^G\/££)^^F^^_^+_g^_(A/£s)_ 
(V^s)8 + 67 iV^sy + • • • + 5o 



[l + 0{s-')) 



where the coefficients bj are non- negative and bo > (but the aj might be zero or negative) . 
Using the bounds ( |5.15| ), one sees that £x'^{s) can be estimated by 



e \xi(s)\ < c 



Vis iiVEsf + 1) 



A scaling argument shows that 



- e Kis) 



ds < c 



Jbi iVesf + 1 Jo 



proving ( |5.16| ). To derive ( |5.17| ), we use that 



r + 1 



< 00 



(5.20) 



(5.21) 



and obtain 



e\xe{s)\ < I £\x',{t)\dt + e\xe{bi)\ 

Jbi 

£ \x (s)\ f'^ \ f r \ 

-^-^Y^ds < — ( / e\x',{t)\dt + £\xe{bi)\\ds 

: ) ( / ^14(01 (^t + ^\xe{bi)\] ds 



d (\ 



bi ds\s , 



According to (l5.2C| ), the inner integral diverges at most logarithmically as s — > cxd. There- 
fore, integrating by parts gives no boundary terms at infinity. 



fci 



i^d. < f 

^ Jb\ 



iKMds + ^\^^^^^\ 
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h (x) 



(h')Jx) 





Figure 2: Typical plots for h'{x) and h'^^{x) in the case cu > m. 



The integral on the right was estimated in (^.211). The last summand is computed to be 



\Xeibi) 



\x'{bi)\ dbi 



de 



and this is bounded uniformly in e in view of ( 5.18| ) and the fact that 61 is smooth in y/e 
and bounded away from zero, (^.71). This completes the proof of Lemma ^?^. ■ 



The above estimates are illustrated in Figure 0, where h' and /i^ are plotted in a typical 
example. The dashed curve describes the asymptotics near x = 0; it is the graph of h'{x), 
where for x{u) one uses the approximate formula x ~ 4:\/]3{y/u — ^/uo) , obtained by setting 



£ in (|5.5| ) equal to zero, dropping the error term and integrating, (4 



In the case uj < m, we fix Umin independent oi uj G {m — 6,m). By choosing Umm large 
and 6 sufficiently small, we can arrange that the function g has exactly one zero on the 
half line (umin; cxd). We set uq equal to this zero. 



giuo) 



. 



Clearly, uq depends on uj. The variable x{u), (4.8), is positive for u > uq and negative on 
the interval {umm,uo)- We set Xmin = 2;(iimm)- The following lemma is the analogue of 
Lemma |5.1| for u < m. The method of proof is also similar, but the pole of /i at x = 
makes the situation a bit more complicated. 

Lemma 5.2 There are constants c,6>0 such that for all uj G {m — 6, m) and x > Xmin, 

< T{x) h'{x) < + -^ (5.22) 

\h"{x)\ < ,^ / . ... + ^ (5.23) 



l/i' (x)| dx < 



yX Xjyiin ~\~ ij 
c 



m — UJ 



(5.24) 



Furthermore, for every xi > the constants c, 5 > can he chosen such that for all 

UJ £ (m — 6,m) and x £ [— xi , xi] , 



h{x) 



1 1 I 2 I 
- log |e x| 



< c 



(5.25) 
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Proof. We now set e = m^ — w^ and p = 2y|7^. Then 



pHu) 



u 



with /3 > 0. Since p{uq) = 0, 



no = ^(l + 0(e)) 



and furthermore, the functions h and p have the expansions 

e |u — Wo' 



/i(u) 

/9('U) 



2 ^°^V(w + m)2 



l + 0(u-i)) 






(5.26) 

(5.27) 
(5.28) 



Since global bounds for x{u) and u{x) would be more difficult to obtain than those in 



Lemma 5.1, we here construct the bounds piecewise. We set Au = iio/2. By decreasing 
5, we can arrange that no/4 > Umiri! and furthermore we can also make the error terms 



in (5.27) and (|5.28| ) as small as we like. Thus we may assume that 

< p{u) < 3 . / — -y/lu — uol for I u — uo I < Ati. 

V ^0 



n \ \/\u-Uo\ 



Integrating from uq to u gives 

£ , I 3 I / M / £ 3 

u — uo\2 < \x{u)\ < 2. — n — no 2 

Uo ■ V ^0 



for \u — uo\ < Au. 



Taking the inverses, we obtain for |u — uol the bounds 

v{x) < \u{x) — uo\ < v{x) for |x| < Ax, 

where we set 



v{x) 
Ax 



x^ Uo ^ ^ 



4£ 

£ 



v[x 



x^ Uo ^ ^ 



(Au)t = 2-^^^/euo = 2-^ (3 £-^2 {I + 0{£)) . 



(5.29) 

(5.30) 

(5.31) 

(5.32) 
(5.33) 



If X > Ax, the second inequality in ( ^.30 ) shows that u — uo > 2 suq, and thus in this 
region there are constants oi, 02 > such that 



oi \/e < p{u) < a2\f£ . 



(5.34) 



Hence for x > Ax, 



1 



du 
< — < 



02 \/e dx ai \fe 
Integration shows that, possibly after decreasing a\ and increasing 02, 



u(x) < u(x) < u(x) 



for X > Ax 



(5.35) 
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with 



u{x) 



02 \/e ' 



u{x) 



ai ^/e 



(5.36) 



If on the other hand x < —Ax, we see from ( 5.3C ) that u — uq < —2 s-uq, and thus p can 
be estimated by 



— u 2 < p{u) < —u 2 
with fei, 62 > 0. We integrate from Umin to u, 

and solve for u. This gives 

u{x) < u{x) < u{x) for Xmin < X < —Ax 

with 



(5.37) 



u{x) 



+ V Wmir 



u(x) 



+ \/^i^ 



(5.38) 



(5.39) 



02 / V Ol 

For any xi > 0, we can, by choosing 5 small enough, arrange that Ax, ( 5.33| ), is 
greater than xi. Thus u{x) is on the interval [— xi,xi] bounded by ( 5.31| ). Substituting 
these bounds into ( p. 271 ) and using ( |5.26| ) gives ( 5.25| ). 

To show that T{x)h'{x) > 0, note that 

du 



h'{x) = h\u) 



dx 



where ^ = p~'^ is positive, and the sign of h'{u) is obtained from ( 5.27 ). 

For the derivation of the inequalities ( ^.22D -( ^.24 ), we consider the three regions x < 
—Ax, |x| < Ax, and x > Ax separately. We begin with the case |x| < Ax. For h'{x) and 
h"{x), we have again the bounds ( 5.11| ) and (5.12), respectively. Using that u > uq/2 as 
well as ( 5. 29] ), ( p.31| ), and (5.32), we obtain 

\h'{x)\ < -^^ < , ,-- . 



v^p^ 



C 3 

3- \u-Uo\~2 < 



\h"ix)\ < 



u^p^ 



+ 



V? p"^ 



< 



Un e^ x^ 



+ 



UqE X 

c 



^0 



§ 4 4 
Un £3 X3 



and in view of ( ^.26 ) and ( |5.33| ), this implies ( ^.22 ) and ( 5.23| ). To compute h^{x), we 
again use ( ^.13|) , but now with 



s{u) = . — |u — uol 2 r(n — uo) 
V ^0 



for \u — uo\ < An , 



where as before F is the step function F(t) = sgn(T). The first summand in ( |5.13| ) is 
computed as follows. 



u{s) 

h{s) 
ehe{s) 



— - — 
Uq +U^ e 3 S3 T{s) 



for Isl < Ax 



(5^6) l(^+(^^,2)ir(s))(l + 0(6)) 



(5.40) 
(5.41) 



(5.27) 1 



log 



14 2 
P3 £3 S3 



2 V/3 + (/?es2)tr(s) 
[2 1 {(3es^)lT{s) 



^3 G p + {pes^)3 r{s) 



l + eO(s°) + e3 0(s3 



l + eO(s°)+e3 OisI 



log{uj + m) 
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Differentiating with respect to s gives the bound 

e|^e(s)| < ces s 3 +ce3 53 , 

and thus 

/■^^ 1 2 2 

/ e\h'^{s)\ds < ce3 (Ax)3 +ce3 (Ax) 

Using ( 5.33| ), we conclude that the total variation of £h£{s) is bounded uniformly in e. In 



order to estimate the total variation of the second summand in ( ^.131 ), we first compute 
x'(s), 

x'{sf = {p{s)u'{s)f = l^^(l + eO(.0) + elo(.i) 

9 u[s) \ 



This is uniformly bounded from above and below, proving ( 5.15| ). Differentiating with 
respect to e using ( ^.41 ) gives the estimate 



e ke(s)| < ce3 S3 . 
Since x{s = 0) = for all e (from ( 5.40D ), integration yields that 



e|a;£(±Ax)| < 
^-e|4(s)| 



±Ax 



e \x'i,[s)\ ds 



1 5 (5-33) 

< cE'i (Ax) 3 < ce 



ds < c£3 (Ax) 3 < c 



Ax 



Furthermore, 



Ax 



£ \Xe{s) 



ds < 



Ax 



Ax 



Ax S \Jo 



e|x^(t)| dt] ds. 



(5.42) 

(5.43) 
(5.44) 

(5.45) 



Using ( 5.42| ), the inner integral is for small s bounded by a constant times ts. Thus when 
integrating by parts, we get no boundary terms at s = and obtain 



Ax 



e \xe{s) 



ds 



Ax 



1 rAx !• 



Ax 
Ax ^ 



ds < c 



(5.46) 



where in the last step we used ( 5.42| ) and ( 5.44| ). Combining ( p. 44 ) and ( 5.46| ) with the 
estimates ( |5.22| ) and ( |5.23| ) , we conclude that the total variation of the second summand 
in ( p.l3| ) is bounded uniformly in e. This shows that ( p. 24 ) holds if the integration domain 
is restricted to x € (—Ax, Ax). 

In the case x > Ax, ( |5.34| ), ( |5.35[ ), and ( |5.3£| ) yield, again using ( |5.11| ) and ( p.l2[) , 



\h'(x)\ < 



c c 

< 



v?p'^ 



1 



M c 

x^ ->/e ~ X a/e Ax ~ x 



c 
< - 



U^p° U-^p^ X* £ X3 ^y£ X'^ 



proving ( |5.22| ) and ( ^.23| ). To compute the total variation of /ie(x), we apply ( 5.13| ) with 



s{u) 



3^/2 



for u > uq + Au. 
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Using that 



we see that 



s{uq + Atij = s( I ^-^^ Ax 



u{s) 



3^/2 



for s > Ax. 



Moreover, from ( ^.27 ) and ( 5.2g| ), 
h{s) 



1, / e 3\/2s-Pe-^ 
t: log 



eheis) = 

e\K{s)\ < 

e\h'^{s)\ds < 



(a; + m)2 3^5 ^ v y v v y; 



1 i + ^^fli 

2 \ 6^/2s-2/3e-5 



r (l + eO(/) + V^O(5-i)) 



/es^ 



< c 



x'(s)2 = (p(.)n'(s))2 



72 1 



/3 



3\/2s Ve 



(l + eO(s") + V^O(s-i)) 



Since s > Ax, we conclude from ( |5.33D that ( 5.15| ) holds. Differentiating the last relation 
with respect to e and integrating gives 



e 14(5)1 < 



l^iMd, < 



s v^ 
c 



< c 



Ax ^ 

00 



e \x^[s) 



Ax 



^/e Ax 
ds < / ^(/ e|3;'(OI t?* + e3;£(Ax) ) ds 



(5.47) 



Ax 



Ax 



where in the last step we used ( 5.43 ) and (|5.47] ). This proves ( ^.24 ) if x > Ax. 
Finally, if x < —Ax, the bounds ( 5.38| ) and ( ^.39 ) give 



(5-11) c 8c 1 

C 



< 



X XYain ~r 02 y "Uinir 



^2 \At 



< 



|/i (x)| ^ ^r^ + ^^ < - < 



C 



This concludes the proof of ( 5.22| ) and ( |5.23| ). In order to prove ( ^.24] ), we apply ( 5.13| ) 
with 

•5(^) = o yP^ ~ 2^^ ^"-"^ '"min ^U <Uq — Au. 

Similar as in the case x > Ax, 



u{s) 



P 



(S + 2Ax)^ for Smin = h V(3 Umin " 2Ax < s < -Ax (5.48) 
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{h')<.{x) 



^rain 



^min 





Figure 3: Typical plots for h'{x) and h!^{x) in the case u < m. 



his) 
ehe{s) 



log 



1 



^2 (1 + 0{e) 



e 4(s + 2Ax)2+ /320(e) 



e") (l + 0((s + 2Ax)-2))') 

eO((s + 2Ax)-2) 



2 /32(i + 0(e))-4e(s + 2Ax)2 



+ eO(s") + 



2 /32 - 4e (s + 2Ax)2 
e|/ie(s)| < ce(s + 2Ax). 



l + 0((s + 2A2;)-2) 
+ eO(s°)+eO((s + 2Ax)-2) 



Integration yields that 

j- — l\X 



(M3) / 2 1 

£\h'i,{s)\ds < Ce I Ax — -/?Umin 



(|5.33|),(|5.2eL 

< ; c e (no - Au - Umin) < c 



Moreover, 



x'is? 



{p{s)u'{s)f = u{l-^e{s + 2^x)A{l + eO{s^) + 0{{s + 2^x)-^)) . 



Using that {s + 2Ax)^ < (Ax)^, we conclude that x'{s) is uniformly bounded from above 
and below. We differentiate with respect to e and integrate to finally obtain similar 
to (|5l^)-(CT), 



e|4l < ce(s + 2Ax)^ 



-^^ e \x'^{s) 



mm 

Ax 



s + 2Ax 



e \Xp[s) 



ds < Ce (uo — Au — Uniin) < c 



(5.49) 
(5.50) 



ds < 



(s + 2Ax)2 

e Ixf (— Ax 



'Ax 



1 



(s + 2Ax)2 V7-A. 



e \x'^{t)\ dt + e |xe(-Ax)| ds 



< 



Ax 



+ 



1 



+ 2Ax 



-Ax 



e|x^(t)| dt + 



""^l-^(»'l<i. <c 



s + 2Ax 



where in the last line we integrated by parts and used ( 5.43| ) and ( ^.49 ), ( 5^ ). This 
yields (5.24) and completes the proof of Lemma 5.2. ■ 



The above estimates are illustrated in Figure y. The dashed curve is the graph of (3x)~^; 
it is the x-derivative of the asymptotic function ^ log(e^x) which appears in ( 5.25] ). 
The next lemma controls the behavior of Xmin- 



25 



Lemma 5.3 There are constants c,6 > such that for all to S [m — 5, m), 



Xmin + 4.l3{m^-UJ^) 2 

Proof. According to the definition of x and ( 5.2^) , 



< c 



< 



m — id 



-2 / p(u) du 



/MO 
--■mill 



g .^ ^Ql (l + 0(n-i))rf^ 



u 



(5.51) 
(5.52) 

(5.53) 



and a calculation using ( 5.26| ) shows that the leading contribution in ^/e to this integral 
is 4/?e~2 (this can be readily verified using Mathematica) . This proves ( ^.511) . Differenti- 
ating ( p. 531 ) with respect to lo and estimating the resulting integral gives ( |5.52| ). ■ 



6 The Region uo > m 

We turn now to the planar equation ( f4.12| ). Consider the two solutions {-d^ ' ,L^ >), b = 1,2, 
with 

lim (i?W,lW) = (0,0) , lim (t?^^),^^^)) = (^,0) (6.1) 



x^—oo 



x— >— oo 



and define in analogy to the transmission coefficients, the quantities (i?oo , -^oo ) by 

\L^£) = lim(^('')(x),L(^)(x)) . (6.2) 



The next lemma expresses the coefficients tab in the integral representation (1.2) in terms 
of these "transmission coefficients." 



Lemma 6.1 The coefficients tab, ( 3-i5 ), are for lo > m given by 

hi = t22 



tl2 



t21 



1 

2 

1 

- e 
2 



-i/3o 



tanhz , 



where 



Po 



arctan ■ 



mri 



r(l)_ r{2) 



+ i 



^^^ - t?S^ + vr 



(6.3) 
(6.4) 

(6.5) 
(6.6) 



Proof. According to ( 4.1C| ), ( fl.5D , and (p]), the functions {'&,L) correspond to the 
2-spinor 



(i/3 
e"^ 
e-2 



h 

e 2 cos ■ 



-6 2 cos 



_h 

i e 2 sin 



x+-d 



x+'d 



I e 2 sm 



x+'d 



(6.7) 



In the limit u -^ — cxd, the function /i ^ (cf. ( |4.9| ) and observe that lim^^n fii") = oj = 
\mir\r^ g{r) according to (|j) ^^^ (^^),( ^34|) ), and x coincides asymptotically with Q.qu, 
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up to an (irrelevant) additive constant. Thus comparing (qJ) with ( p.?]) and using (|6^ 
gives 



e 2 \ 



f(i) _ f e 2 ] f(2) _ 

\ — e 2 / \ — ^e 2 



. _1^ 
ze 2 



Hence the fundamental solutions Xq and Xi, which are characterized by (|2.q) , are the 
linear combinations 



X, = ie^(X«+^x(2)) 



X2 = ie-^(-X«+.x(^)). 



We next consider (|6J) in the limit u -^ +00. According to (|49|), (|j), (|j)-(gj), 
and (|1|), 

1 7 / \ 1 , w — m 

lim ft(x) = - log = — zB . 

x^oo 2 a; + 7TT- 

Also, /3 goes to zero in this limit. Hence using (2.12) and (6.2), one sees that 



/ 



, 1 ±iio 

/oo 1/2 — TJ e 2 



V 



ibexp 
ibexp 



r(l) , ,-„q(1)' 



r(l) ,-,q(1)' 



+ i exp 
+ i exp 



-t^co + i't'oo \ 



r(2) •n(2)^ 



Substituting this last formula into ( 2.14 ) yields 



i/30 



*i/2(a) = e^ 2 



le 2 sm a + 






±e 2 sm a + 






A short calculation shows that 

|ti|' = \t2?. (6.8) 

Together with (2.15), this immediately yields ( |6.3D . Furthermore, it is obvious from ( |2.15| ) 
that ti2 = ^21- Thus it remains to compute ti2. According to ( 2.15 ) and (|6.8|) , we have 



tl2 



1 1'^^ tit2 



2tt 







2|i2p 



1 /"^^ h 



Att 



t2 



-if3o /•27r 



ipi sin(a + <y9i) + p2 sin(a + 992) 
4vr Jq i/9i sin(a + v^i) - /32 sin(a + </72) 



da , 



where we set pi = L)^/2 and ipi = t?S/2. It is convenient to shift the integration variable 
by a — > a — (/92 and to divide the numerator and denominator by p2- This gives 



-i/3 [■2tt 



tl2 



ip sin(a + 93) + sin a 
4vr ^0 ip sin(a + 9?) — sina 



(ia 



with p = p\l P2 and 93 = (^1 — (^2- We express the trigonometric functions as exponentials 
and set p = pe~^'^, 



tl2 



i/3o /■2^ (^ _ ^) g2m -(^-i) 



4vr 7o (/" + i) e^*° - (/^ + 



da . 
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Setting z = e^*", the a-integral can be regarded as an integral along the unit circle in the 
complex plane; more precisely, 



-i/3o 



tl2 



{^ — i) z — {/J, — i) dz 
47ri 7|^|=i {jl + i)z-{iJ: + i) z 



(6.9) 



This contour integral can be computed with residues as follows. According to ( |6.1| ), 
limj;^_oo "doc — t?oo = — vr. A comparison argument using the differential equation for 



1?, ( 4.12| ), shows that i!}^^> — -d^"^' takes values in the interval (— 27r,0) for all x. Hence 
— TT < 93 < 0, or equivalently, 

Im fi > . 



As a consequence, the integrand in ( |6.9D has only one pole in the unit circle, at z = 0. We 
conclude that 



2 fi + t 



1 

- e 
2 



-i/9o 



>/i - y/i/fx 



and this coincides with (|6.4|) . 

The following two lemmas control the behavior of {'d, L) for large x. 

Lemma 6.2 There is c > such that for all uj £ {m, m + 5) and x £ (0, cxd], 

|t9(x)-??(0)| < c, \L{x)-L{0)\ < c. 

Proof. According to (^j|), there is xq > such that 

1 + /i' sin(x + "i?) > — for all x > xq. 



On the interval [0, xq], we can control "i? by integrating the iJ-equation in ( [4.12 ) 



|i?(a;o)-i9(0)| = / h' sm{x + ^) dx 
Jo 

In the region x > xq, we again integrate the equation, 

PX j-X 

^{x) - ^(xo) = / h' sin(x + ^)dx = - - 

Jxu Jxn 1 + h' 



(5.1) 



/l' 



CXo 



d 



(6.10) 



(6.11) 



(6.12) 



sin(r + -d) dr 



(cos(r + ■&)) dr . 



We integrate by parts and, using ( 6.11 ) and Lemma ^^, we find 

|i?(x) - i?(xo)| < 2 (|/i'(2;)| + \h\xo)\) + 4 r (|/i"| + /i'^) dr < c . (6.13) 

The second statement in ( |6.10 ) follows similarly by integrating the L-equation in ( [4.12| ). ■ 



Lemma 6.3 There is c > such that for all uj S (m, m + 5) and x G (0, 00], 



\Mx)-Mo)\ < 



UJ — m 



\L^{x) - L^m < 



UJ — m 



(6.14) 
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Proof. Differentiating through the ODEs in ( [4.12| ) with respect to uj gives 



^', 



L' 



h' cos(2; + '&)^uj + h'^ sin(x + ■!?) 

L' ^^ + hi sin(x + ^) 

— h' sin(x + t?) i?cj + /i^ cos(2; + -&) 



(6.15) 
(6.16) 



The differential equation ( 6.15| ) can be solved using the method of variation of constants. 
The solution is 



^u.{x) - ^^(0) = e^(") r e-^(^) h'^{T) sin(r + {}) dr . 

Jo 

Lemma |6.2| yields that 



(6.17) 



Jo 



and the estimate ( ^.3| ) in Lemma ^[^ gives the first part of ( |6.14 ). To derive the second 
part, we integrate ( |6.16| ) and apply again the integration-by-parts technique of Lemma |6.2| , 



h'^.. 



-— cos(r + -d)] dr + I h'^ cos(r -|- ■&) dr 
d-r J Jo 



\^^(-)-^^m^ J^ l + /.'sin(r + ^) 

< 2\h'^^\{x) + 2\h'^^m + 2 [ {\h" ^^\ + h'^ \^^\ + \h' ^'J) + [ \hl\. (6.18) 

Jo Jo 



Using the estimates of Lemma |5.1| and Lemma 6.2, the only problematic term is the 
integral of | /i' i?^ | . But from ( 6.15| ) and ( [4.12| ) we have 



\h'^'. 



\h" sin(x + 1?) T?^ + h' h'^ sin(x + d)\ 

1 c 



•12 



< \h"^J + \h'h'\ < 



UJ — m (1 + x)^ 

where in the last step we used ( |5.1| ) and the first part of ( |6.14| ). 
We remark that by combining ( |6.16| ) and (4.12), we can write the L^-equation as 

L'^ = -^'^u + hi cos(x + ^) . 



Although this looks very similar to (6.15), it seems difficult to deduce the second part 
of ( |6.14| ) by integration (note that the total variation of i? need not be bounded uniformly 
in io). This is the reason why we instead used an integration-by-parts argument. 

We are now in the position to prove that in the integral representation (|1.3|), the 

5 

contributions for a; > m decay in t at least at the rate t~ 6~^. Consider the two fundamental 
solutions (i?(^), L^^-*), ( |6.lD . For negative x, the function h'{x) is smooth in co. Furthermore, 
h'{x) is computed to be 



^'W = '='('•) ^ ^ '=''.'('■) ,._, 



du 



du dx 



r^ + a? dx 



(6.19) 



Using that A decays exponentially as n ^ — oo, and that for large negative x, u'{x) is 
bounded away from zero, we see that h'[x) decays rapidly as x — > — oo, locally uniformly 



in oj. Thus standard Gronwall estimates applied to the differential equations ( [4. 121) yield 



29 



that (^'^'(0), L(^-'(0)) depends smoothly on uj. Hence Lemma S^ and Lemma 6^ give us 
information on the transmission coefficients, namely 

|^W|,|^W|<c and |a^^W|,|a^^W| < ^_. (6.20) 

uj — m 

Next we consider the propagator ( |1.3D for x in a compact set K and ^o with compact 
support. Again, standard Gronwall estimates starting from the event horizon yield that 
the fundamental solutions ^^"^"(x) depend smoothly on u), uniformly for x ^ K. Hence 
the only non-smooth terms are the coefficients t^^"'- According to Lemma p.2| , these 
coefficients have the same regularity as the transmission coefficients, ( |6.20 ). Furthermore, 



Theorem 2.1 allows us to restrict attention to a neighborhood oi oj = m, and thus we 
may assume that the square bracket in ( |l.3| ) has compact support. We conclude that this 
square bracket satisfies the assumption of Lemma O] (with a = a; — m), and thus its 



Fourier transform decays like t e" 



7 The Region uj < m 



For ui < m, the coefficients tab in the integral representation (|l.3| ) have the simple explicit 
form ( 2.15| ), and thus our task is to analyze the w-dependence of ^j^'^"(x). We again work 



in the variables {'&,L) and set 

(j){x) = X + t?(x) . 

Recall that ^j^'^"- is the fundamental solution with exponential decay at infinity. The 
following lemma shows that this implies that lima;^oo (pi^) = — oo. 

Lemma 7.1 There is a constant C independent of lv such that for all x > 0, 

<t>{x) < C -logx . (7.1) 

Proof. Using the bounds ( p.22| ) and that h' is positive, we have 

^' > l-^e-'f' (7.2) 

with c independent of u. Suppose that ([7.ip were false for some x = xq and C = logc. 



Then (p\% ) implies that 

c/>'(x) > I . (7.3) 



Hence aX x, (p is monotone increasing, whereas the right side of ( |7. 1| ) is monotone de- 
creasing. As a consequence ( [7.1[) is violated on an open interval (x,x + e). Furthermore, 
by continuity ( |7.1| ) is violated on a closed set. We conclude that ( |7. 1| ) is violated for all 
X G [xo,oo). This means that ( |7.3| ) holds for all x > xq, and integration yields that 

lim (j){x) = oo . ('''•4) 

x— >oo 



To finish the proof, we shall show that (|7.4D implies that the corresponding two-spinor -0, 
( [4.13| ), grows exponentially at infinity, giving the desired contradiction (note that since 



i^kujn (jgga^ys Q^^ infinity, X, X, and tp also vanish at infinity, see (|2.2|), ([4.1D and (7.17)). 



According to (7.4), ip behaves for large x asymptotically as 

V; = e^ I "1^. 1 (1 + 0{e-^'^)) 
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Furthermore, using ( [4.12 ), 



{(p-L)' = 1 + h' (sinh (/> - cosh (/>) = 1 - h' e'^ " 1 + 0{e-^) . 
Hence for large x, (j) — L ^ x, and so ^ grows asymptoticahy hke ^ ~ e^. ■ 

The inequality (^]^) shows in particular that 

(/)(x) < — 7T for all X > xi (7-5) 

with xi = exp(C + ^). We next show that (/) leaves the region {(f) < — i} for positive x. 
Lemma 7.2 There is xq > v' > with v independent of uj such that 

4>{xo) > -\ . (7.6) 

Proof. We introduce for x > the function 

1 1"^ 
l{x) = log^ / h'{T)e~^ dT . 



Since by (|5.22| ) the integrand is positive, t9 is monotone decreasing. According to ( 5.22| ) 
and ( |535| ), 



lim ^{x) = — oo , lim ;d{x) = oo 

X-+00 X— >o 

and so there is a unique xq with 

Now, choosing < y < z, we have 

/i'(r) e^^ dT > e-^ {h{z) - h{y)) . 



Using (|5.25|) , we see that for small y, 

/•oo 

/ h'{T)e-^dT > 4, 

implying that xq is bounded away from zero, uniformly in cv. 
We shall now prove that i?(x) is a lower bound for t?, i.e. 

??(x) > ??(x) for all X > xq. (7-7) 

Thus in the region x > xi, we apply (|7.5D to get the estimate 



1 ./ „-<* 1 ./ 

= n sum (p ^ — ^ -- ^ — 

We separate variables. 



t9' = h' sinhc/. < — /t'e"-^ = — /i' e-(^+'') . 
4 4 



^i9\'/„\ ^ ^ i,// 



(e'^)'(x) < --/^'(x)e-^ 
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and integrating (using that e^^°°' =0), we find 



Mo=) > 



1 f°° 

-j^ h'{T)e-^dT. 



Thus •§_ is indeed a lower bound in the region x > xi. 

It remains to show that i9 > ?9 on the interval [xo,a;i]. Let us assume the contrary. 
Then {} and '& meet somewhere on this interval. Let 

y = sup {x I ■i?(x) = T?(x)} . 
Then ^{y) = ^{y) < -|, and thus 

This contradicts the fact that i!}{x) > ??(x) for all x > y. ■ 



The next lemma controls the behavior of (j) near the origin and "matches" the solution 
across the singularity at x = 0. 

Lemma 7.3 Suppose that for given ki < and K2 > 0, 

-- < 4>{x) < for allx e {ki,K2). (7.8) 

Then there is ki with ki < 0, ki < ni and a parameter A > such that 

C _A e'^^^) + r(x) < ^(x) < -Ae^(^) for ki < x < 

\ -Ae'^^^) < (/.(x) < -Ae'^(^)+r(x) forO<x<K2 

with 

r{x) = e'^(^) r e-^^^^ dr . (7.10) 



JO 

Note that the function r(x), ( [7.101 ) is finite according to ( ^.25 ) 



Proof of Lemma |7.5| . Let us first assume that Ki < 0. We set ki = ki. We choose r] 
with < r] < min(— Ki,K2). For negative x, (j) satisfies according to ( [4.12| ), the equation 
4>' = 1 + h' sin (f). Using the bounds (|7.8| ) as well as the fact that h' is negative, we obtain 
that 

h'{ct> + ^^) < (/)' < l + h'ct> (7.11) 

for all X S (ki,0). We choose x in the interval {ki, —rf) and consider the inequality (|7.1l| ) 
on the interval (x, —rf). The inequality on the Ihs can be solved by separation of variables 
and the rhs by variation of constants. This gives the explicit bounds 

gh(x)-M-'y)0(_^) + gM^) / e-'^Wdr < 0(x) < -^^ for -k < x < -r/ (7.12) 

.hr I- a 



with 



a = e'^(-)-M-'7) ^t;}) . (7.13) 

l + (/>(-T?) 
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If X is positive, then according to ( 4.14 ), (p satisfies the equation (f)' = 1 + h' sinhc/). 
Using (|7.8|) and that h' is now positive, we get the bounds 



l-h'cl) < -(/)' < -h' {(j) - (j)^) 



on (0,x). 



(7.14) 



We choose x in the interval (??, k) and integrate these bounds from r] to x. This gives the 
bounds 



/3 



1 + /3 



Jr] 



for 77 < X < K (7-15) 



with 



/3 



^h(x)~h(vi) 



<Piv) 



We now show that 



1 - <^(r?) 

= \im (j){r]) . (7.16) 

Consider </)(— r/). From ( 7.12| ) and ( [7.8[ ) we have for fixed x in the interval —k < x < —rj, 

< 0, (7.17) 



lira Si-rj) 



1 

- < 
2 - 



a 



(x) < 



1-a 



and thus there is some ao < for which a > ao if ?? is sufficiently small. According 
to (|]2|), the factor e''(^)-'^(-'?) in (|7.13D tends to +00 as 77 \ 0. We conclude from (|7l|) 
that 



lim 



0, 



r,\0 1 + (t>{-v) 

implying the Ihs of ( [7.16| ). A similar argument using the rhs of ( [7.15 ) gives the rhs of ( |7.16| ). 
Since the planar equation ( [4.6| ) has smooth coefficients, it is obvious that ip{u) is 
smooth and non-zero. Using the ansatz' ( [4.10D and ( 4.13 ) as well as ( |7.16| ), we see that 
the following limits exist, 



L(-,))-fc(-7,) 

lim e 2 

'?\o 



M-v) 



. L{r])-h(r]) 

H-V) = V'U'=o = lime 2 

ri\0 



-Hv) 



<t>{v) 



We consider the two cases lim^\o(-^(^) — ^iv)) = ^^'^ ¥" separately. In the first 
case, the second components must have a non-zero limit (because V'(O) ¥" 0)) ^^^ thus 
lim^^o e~'* </'(=t^) = — cxD. In the second case, the limits lim,y\^o e~'* V(=t'?) must 
exist and be equal. We conclude that 



lim e-'^^-^) (/.(-r/) = -A = lim e'^^^^ 0(ry) 

»?-+0 r;— >0 



(7.18) 



for some A G [0,oo]. In the case ki = 0, this matching of the two ansatz' shows that 
(p < for negative x, and thus we can make ki slightly negative and repeat the above 
construction. Again using (|7.17| ), one deduces that A must in fact be finite. We finally 
take the limit rj ^ in ( [7l^ ) and ( [tIsD to obtain (JTj). ■ 



Our next goal is to bound i9(xinin) uniformly in uj. To this end, we combine the a-priori 
estimates for large x (Lemmas [7.1| and Lemma |7.2| ) with the estimates in a neighborhood 
of X = (Lemma [7.3| ). For negative x outside of this neighborhood we can use similar 
methods as in Lemma p^. 
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Lemma 7.4 There is c > such that for all uo G [m — 5, m) 

\-&{Xmin)\ < C . 



Proof. First let us verify that the assumptions of Lemma 7.2 are satisfied for a particular 



choice of ki and K2- To this end, observe that (l){x) has no zero for x > 0, because otherwise 

(p'ix) = l + h'{x) smhcj) = 1 , 



violating the fact that (/>'(x) < at the largest zero (recall that Lemma |7. 1| implies that (j) 
is negative for large x). Thus 

(pix) < for all X > 0. (7.19) 

As a consequence, sinhc/) < 0, and thus using ( ^.22 ), 



b'{x) < 1 + h'{x) (/) for all x > 0. 



Integrating this inequality from a given positive x < xq to xq and using (7^), we obtain 
the lower bound 

1 r^o 

(t>{x) > l{x) = — e'^(^')-'*(^'') - e'*(^) / e-'^W dr for < x < xo (7.20) 



(this is indeed quite similar to the second part of ( 7.15| ), but now we have solved for 



at the lower limit of the integration range). According to (5.25), linix^Q (l)(x) = 0. We 
conclude that the assumptions (|7.8|) are satisfied for ki = and K2 > sufficiently small. 
We can further decrease ki and increase K2, provided that the bounds in ( [7. 9] ) all take 
values in the strip (— i, 0). 

The parameter A in ( |7.9| ) can be bounded a-priori. Namely, were A sufficiently large, we 
would get a contradiction to (|7.20| ), whereas a very small value of A would be inconsistent 



with (^^q). Thus we can find parameters < Amin < Amax such that 

'^min \ A ^ Ama.x . 



min 



As a consequence, in (|7.9|) the lower bound for A = Amax and the upper bound for A = A 
are a-priori bounds for (p. We choose X2 such that these bounds take values in the strip 
(—^,0) on the interval [a;2,0). Then we have a-priori bounds for (j){x2), and thus also for 

^{X2) = </'(x2) -X2, 

^min < ^(X2) < l9max • (7.21) 



The bounds ( [7.21 ) are uniform in u). This is not surprising since the differential 



equation for •& involves only h', which according to Lemma 5.2 is bounded uniformly in 
uj. To see this rigorously, one must be careful because Amin and A^ax do depend on lo. 
Namely, according to ( |5.25| ), h involves the additive constant ologe^, which diverges as 
ui y m. this implies, according to (|7.18| ), that 

2 2 

can be chosen uniformly in e. Using these scalings in ([7.S|), one sees that the estimates for 
1? and X2 are indeed uniform. 
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It remains to control "d on the interval [xmin,X2]. According to ( ^.22 ), there is i? > 
independent of uj such that 



for X e [Xmin + R,X2 — B\ 



1 + h'{x) sin(x + '&)>- 

(note that this last interval is non-empty in view of Lemma [5^ ) . On the bounded intervals 
[a^minj a^min + R) and (x2 — R, X2\ wc Can control i? directly by integrating the equations in 



a method similar to ( 6.12 ). In the intermediate region, we integrate by parts and obtain 
similar to ( |6.13| ), 



\^{Xm^n + R)- ^{X2 - R)\ < 2{\h' {x^^ + R)\ + \h' {x2 - R)\) + 4 {\h"\ + h'^) dr , 



X2-R 



and the terms on the right are all uniformly bounded according to Lemma |5.2| . ■ 

The next lemma controls the w-dependence of 'Q. 

Lemma 7.5 There is c > such that for all uj £ [m — 5, m), 

m — IjO 

Proof. In the proof of Lemma |7.4| , we have verified that the hypothesis of Lemma 7.3 
are satisfied, and thus '!?(0) = for all oj. Hence "(^tjCa^min) is obtained by integrating the 
differential equation ( |6.15| ) from Xmin to zero. This gives in analogy to ( 6.17 ), 



VujKXTainj 



^(^™) / e-^(^) h'^{T) sin(/)(r) dr . 



By definition of ^5''^", lim„^_oo iv(u) = 1 (see [|, eqn (3.31)] and (|4To| )). Standard 
Gronwall estimates on the interval (— oo, nmin) show that L{x^^^) is bounded uniformly in 
UJ. Furthermore, it was shown in Lemma [5^ that {m — uj) h^ has bounded total variation. 
Thus to finish the proof, it suffices to show that there is c independent of uj such that 



e-^M smdjir] 



< c 



for all r G \x^ 



,0). 



(7.22) 



The integration-by-parts technique of Lemma ^^ yields that L is uniformly bounded in 
the region [xmirn 3^2] with X2 as in the proof of Lemma W^ (for more details see the last 



paragraph of Lemma 7A, where this method is used to estimate t?). On the interval (x2, 0), 
the a-priori bounds for cj), (|7.9|), show that 

|sin</.(r)| < |0(t)| < c e'*^")-'^^^^) (7.23) 

(with c independent of a;). Furthermore, 

(h-L)' = h'{l-cos(t>) < \h'\(f>^ < ca;~5 , 

where in the last step we used ( 7.231) , ( |5.25|) , and ( 5.22 ). Since x^s is integrable, 

ih-L)\l^ < c. 

We exponentiate and use that L{x2) is bounded to obtain 

^-L{r) ^h(r)-h{x2) < g^ (7_24) 

The inequality ( 7.22 ) follows by combining ( f7.23 ) and ( 7.24 ). ■ 
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8 Proof of the Decay Rates 



We are now ready to finish the proof of Theorem 1.1. In view of Theorem 2.1 and the 



consideration in the last paragraph of Section |^, it remains to show that the contribution 
to the propagator ( pT^ ) for uj £ (m — 6,m) has the decay ( |1.6[ ). Since the coefficients tab 
are trivial for u < m, (|2.15| ), the contribution to the propagator ( |1.3| ) simplifies to 



-TT ^-^ ^-^ 



|k|<A:o |n|<no 



-,ku)n( 



a,b=l 



kum 



*n> 



Since ^o has compact support, it suffices to analyze the w-dependence of ^^'^"(n) for u 
in a compact set. 



According to the separation ansatz (2^), we must only analyze the radial function 
X (the angular part Y is clearly smooth in cj). To see the w-dependence of X in detail, 
we substitute ( [4.10 ) into ([4. 51) and (fl.l]) . This gives, exactly as in the case |a;| > m, the 
formula ( |6.7D . For fixed u, the function h in (|6.7D depends smoothly on lo. Using that h 
vanishes at the event horizon (because limr\ri /(^) = i-o = lim^^ri 9(1") according to ( |4.7| ) 
and ( [4.3D , ([4.4| )), our normalization condition for 'I'^'^" near the event horizon ( |2.9D yields 
that 



1 = lim \X{u) 



2 lim e 

U—^ — 00 



-L(u) 



and thus limu^„oo L[u) = log 2, independent of lo. Furthermore, an argument similar 
to ( |6.19| ) shows that h'{u) has exponential decay as n ^ —00. Hence standard Gronwall 
estimates yield that L{u) is bounded and depends smoothly on a;. Furthermore, Gronwall 
estimates in the finite region between Umin and u show that the difference 'd{u) — '(^(nmin) 
is uniformly bounded and smooth in u). Writing 

(X + I?)(n) = (x + ^)(nmm) + ((x(n) - 2;(Ui„in)) + (^(-U) - 1?('Umin))) , 

we conclude that the only possible non-smooth terms in (|6.7D are the factors cos(0min/2) 
and sin((/)mm/2) with 0mm = a^min + ^(a^min)- 

We next consider the factors <,'^f^'^ \ ^o> in ( ^-ll) . Again from Gronwall estimates, 
one sees that for to > m, the expectation values <c'^2'^"' \ ^o> depend smoothly on lo, and 
thus our assumption ( |l.5D implies that 



r2 



lim <*^'^" 

w'\m 



*o> / 



Except for the additional phase factors, the expectation values are smooth even for co < m. 
To compute the phases, we consider ( |6.7[ ) in the asymptotic regime u — > —00, and compare 
with ( p.7| ) and ( |2.8D . This shows that for lo G {m — 6,m), 



-.kum 



^^Ka;n | ^^^ ^ ^^ ^^ g^p 



r2 a2 exp 



with coefficients Oa which depend smoothly on co and are non-zero (indeed lim^y'rn o-a = 
1). Since the factor r2 is non-zero, we conclude that <_'^i^ \ '^q> has a non- vanishing 
contribution which oscillates like exp(z(^mm/2). Using ( ^ ) and ( |6.7| ) in (^?l|), we can write 
the propagator in the region u) £ [m — 5, m) as the Fourier integral 



e "*'^* f si e~*(^™'""'''^(^'""")) + S2 + sa e*^^'™'""'"''^^™'"^M d^i^ 
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with coefficients Sj wfiicli are smooth in u) and S3 7^ 0. According to Lemma p.3| 



Lemma 7.4, and Lemma 7.5, the three contributions to this Fourier integral satisfy the 



hypotheses of Lemma 3^, Lemma 3J., and Lemma |3.3| , respectively (with a = m — u> 



Hence the first two terms decay like t s ^, whereas the last term gives the desired decay 



rate ~ t e. This concludes the proof of Theorem l.L 



9 Probability Estimates 

We now proceed with the proof of Theorem |l.2| . We want to compute the probability p, 
( |1.7D . We begin with the following lemma. 

Lemma 9.1 For any Schwartz function f G S(R x M), let Aj- be defined by 



A, 



t^oo 



—00 J —00 



lim / du duj div' e-'^^-^'^^'^""^ f{uj,Lo') . 



Then 



/oo 
f{LU,L0) doj 
-00 

Proof. We integrate by parts to obtain 

/■oo 

doj I diu' e-'^^-^'^^^^""^ f{io,iu') 

1 



and A^ = 



(9.1) 



00 J —00 

00 /•oo 

duj / duj' 

00 J —oo 

'OO POO 

dio / duj' 



(t±n)2 + l 

1 



(9. + 1)(5.' + 1) e-^(---')(*±")) /(..,a;') 



-00 ./-oo (t±n)2 + l 
where g is the Schwartz function 



(9.2) 



Since the factor ((t it n)^ + 1)"-*^ is integrable in u, we can integrate over u, apply Fubini, 
and use Lebesgue's dominated convergence theorem to take the limit t — > 00 inside the 
integrand, 



A, 



/UQ /'OO /'OO 

du du! did' 

-00 J —00 J -co 

00 /"OO 

dw I duj' g{u>,uj') lim 



1 



(t± 1^)2 + 1 
"0 i 



e-ii^-^')it±^) g(u;,u;') 



t^<^J-oc{t±u)^ + l 



-i{u.~u:'){t±u) ^^ ^9^3) 



In the case "t — u," we introduce a new integration variable a = t — u and get for the 
inner integral 



lim / du- 

t^oo 



(t±'u)2 + l 



-i{w-uj'){t-u) _ -^jj^ 



1 



' t — UQ ^^ I -^ 



-i{w-u:')a^^ = 0. 



This proves that A^ = 0. 

In the case "t + u," we obtain similarly an integral over the real line, which can be 
computed with residues. 



r>«o 

lim / du- 

i— >oo 



1 



-i{LJ — U!'){t+u) 



(t±'u)2 + l 



- — e-'^"^-^')" da = ^e-l'^-^'' 



a2 + l 
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We substitute this formula as well as (|9.2| ) into (|9.3| ) and integrate by parts "backwards" , 

/CO /"CXD 

duj / duj' g{uj,J)T:e-\'^-'^'\ 
-oo J —oo 

du r du'f {uj,uj') ((9^ + 1)(9^, + l)e-l"-'^'l) (9.4) 



vr 

-oo ./ — oo 



A short explicit calculation shows that the derivatives can be computed in the distribu- 
tional sense to be 

(c'^ + l)(c'^' + l)e-l"-"'l = 2 6{oj-oj'). 

Substitution into ( |9.4| ) gives the desired formula for A^. ■ 

We remark that the above lemma cannot be obtained by naively interchanging the orders 
of integration. 

Theorem 9.2 The probability q for the Dirac particle to disappear into the event horizon, 
defined for any e > by 



q = lim / (^7-'^)(t,x) ^■j d/i , (9.5) 

*^°° J{ri<r<ri+e} 

is given by 

q = - Y. E / ^'^ E 4r<<'''' I *o><*o I ^^"> (9.6) 



\k\<ko \n\<no a,b=l 



with 






^ab \ n 4.ku}n j.ku)n „--fl,.l^ ^ ^^■'' 



We remark that p + q = 1, since we know from P] that the probability for the Dirac 
particle to be in any compact set tends to zero as t ^ oo. 



Proof of Theorem 9.i. In the variable u, we need to compute the probability for the 
Dirac particle to be in the region u < uq, where uq may be chosen as small as we like. 
Thus we can work with the asymptotic formulas near the event horizon, with error terms 
which decay exponentially fast as uq ^ — oo. More precisely, a straightforward calculation 
shows that the probability integral in (|9.5| ) coincides asymptotically with the integral of 
the scalar product <.|.> on the transformed spinors (see |Q, eqn (2.15)]). Thus it suffices 
to consider the probability 



du / dcos'd / dip<'^\ '^>\{t,u,i),v) , 

-oo J-l Jo 



(9.8) 



and let t — > oo. Due to the additivity of the probabilities corresponding to the angular 
momentum modes (which are orthogonal with respect to the scalar product <.|.>), it 
sufhces to consider a solution of the Dirac equation with fixed angular momentum quantum 
numbers k and n, i.e. 

-1 />oo 2 

*(t,x) = - / duj e-'^^ y C ^a'^"(x) <^^'^" I ^o> • 

•^ °° a,6=l 
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We substitute this formula for the propagator into ( |9.8D and carry out the angular integrals 
to obtain 



1 f'^O /'OO roo 

[t) = ^ du duj du' e-^^"^-^')* 



" -'— oo ^ — oo J —oo 



f^UQ fOC ^oo 

^— oo 

2 

X 

a,b,c,d=l 



Y, c <^" I ^o> t^; <^'^' I ^o> <^d' I x:>{u) . 



Substituting for X'^ the asymptotic formulas ( p. 51 ), valid near the event horizon, we obtain 
with an exponentially small error term 



1 r^o /"OO /"OO '^ 

^ J —oo J—OO J —oo „ !, „ J 



i{uj—w')u I J^Lj' flil i{u)~LLl')u 



X c <^^ I *o> t^; <^^' I ^o> (/c";/r+ e-*^"-" ^" + /,-:/r- e*^"~" ^" j • (9.9) 

Since we cannot expect the integrand to be smooth when w or lo' is equal to =bm, we must 
use an approximation argument. Namely, the integrand is bounded and has rapid decay 
in u and oj' . Thus we can approximate the integrand in L^ by a Schwartz function, and 



applying Lemma 9.1 we obtain 

^=^/ '^ E '^^ab'^dMfa+ <^t I ^0> <^^ I ^0> . (9.10) 

•'~°° a,b,c,d=l 

It remains to compute the factors /^/tt+. In the case \uj\ > m, we conclude from (|2.8| ) 
that 

Mf:+ = S,i 6ai . (for 1^1 > m). (9.11) 

On the other hand if \uj\ < m, using ( |2.15| ) in ( p.lO| ) we must only compute l/i^P- To this 
end, we again consider ( |6.7D . Using that h vanishes asymptotically near the event horizon, 
one sees that 

lim(|X+|2-|X_|2) = 0, 

r\ri 

and thus l/i^P = |/i^_P- Furthermore, our normalization of the fundamental solutions 
near the event horizon (|2.9|) yields that |/f+P + |/f_P = 1, and thus we conclude that 

l/r+P = I (for 1^1 < m). (9.12) 



Substituting ( 9.11 ) and (|9.12|) into ( [9.101) and using (|2.15|) completes the proof. 



Proof of Theorem ll.Sj . Since the initial data is normalized by <^o|^o> = Ij by taking 
the inner product of ^o with (|1.3|), evaluated at t = 0, we obtain that 



1 /'OO 2 

^ = - Yl E / ^"^ E tar<'^t"''' I 'f OX^O I ^a"''> • (9.13) 



fc|<feo |n|<no ■' °° a,fe=l 



As remarked after the statement of Theorem |9.2| , p = 1 — q. Thus p is obtained by taking 
the difference of ( pT3| ) and (^. Using ( plsD , we get (|r|). 
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For the proofs of (i)-(iv), it again suffices to consider a fixed angular momentum 
mode. Since the energy distribution in the interval [—m,m\ is absent from ( |1.8D , it is 
obvious that (iii) holds. 

To prove (ii), we introduce a vector v^ E C^ by 

< = <vI'^|*o>, a = 1,2 

and remark that in the region |a;| > m we can write the integrands in ( |9.13| ) and ( |f.g| ) as 

{T^ v'^ I v"^) and {A'^ v^ \ v^) , 



respectively, where, using Lemma 6.1 

rpUJ 



1/2 t 



t 



12 



12 

1/2 



and 



A^ 





1/2-2 \t' 



w |2 
12 1 



An easy calculation shows that T > A. Thus from ( p.l3D and (2.15) 
1 /■" . , ,,, , ,,, ,o 1 

TT 
1 

vr 



1 



> - 



duj \<^f I ^^>|^ + - / duj {T^ v'^ I v'^) 

m " JlR\[—m,m] 

m 1 f 

dL<j|<^^|*^>|2 + - / du {A"^ V'^ \ V'^) 

m 



1 r™ 

TT J-m 



(9.14) 

(9.15) 
(9.16) 



and this is strictly larger than p because in case (ii) the first summand is positive. 

To prove (i), we note that the factor |<^2 I ^o>P is positive on a set of positive 
measure (by continuity in w). Thus it suffices to show that 



\-'^\tl2? > 



for all a; G M \ [— m, m\. 



Using the explicit formula ( |6.4| ) in Lemma 6.1, this holds iff 



tanhzl < 1 



with z as in (6^). Using ( 4.12[) together with (|6.lD , we see that — 27r < -i? 



(1) 



(9.17) 

i?S^ < 27r 



(by the uniqueness theorem for ODEs). Then from ( |6.6[ ), 



vr vr 

- < argz < - 



It follows that |e^^ — 1| < |e^^ + 1|, giving ( ^.171) . This proves (i). 

Finally, if ( |1.5D holds, then we saw in ( |8.2| ) that <^f |^o> is non-zero for uj G {m—S, m). 
Thus (iv) follows from (i) and (ii). ■ 

Given the fact that he Fourier transform of a C°° function with compact support is an- 
alytic, one might think that <$^'^'^ | \I'o> should be analytic in uj, implying that the 
cases (ii) and (iii) cannot occur. However, it is not at all obvious that the solutions of 
our ODEs should depend analytically on uo. Should this be the case, one could still make 
sense of (ii) and (iii) by slightly weakening the assumptions on the initial data. 

We conclude by describing the class of initial data for which the Dirac particle must 
escape to infinity, with probability one. 
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Corollary 9.3 The probability p is equal to one if and only if the initial data satisfies for 
all k, CO, and n the following conditions, 

r <^^'^"|^o> = if\oj\<m 

Proof. It again suffices to consider a fixed angular momentum mode. In view of ( |9.16| ), 
p = 1 only if 

/m 
d(j|<^^|^^>|2 = 0, (9.18) 

-m 

and thus the energy distribution of the initial data must be supported in the outside the 



interval {—m,m). Furthermore, the inequality in ( 9.15 ) must be replaced by equality, and 
thus 

(5"^ v'^ I v"^) = for all a; G M \ [-m, m], (9.19) 

where the matrix S"^ is defined by 

^ — I TIj- 9 Uc; |2 
^^12 ^l''12l 



The eigenvalues of S"^ are zero and ^ + 2|ti'2p ^ 0- Hence (9.19) implies that v'^ must be 
in the kernel of S"^, i.e. 

< = -2 tf^ v^ . (9.20) 

Conversely if (|9l^ ) and ( |93o| ) hold, it is obvious from (^ll)-(pl6|) that p = 1. ■ 



One can also understand directly why p = 1 for special choices of the initial data. Indeed, 
to obtain such initial data, one can consider the physical situation where a Dirac particle 
at time t = —oo comes in from spatial infinity. If we take as our initial data the corre- 
sponding ^(t, x) at t = and reverse the direction of time, the solution to this Cauchy 
problem will clearly escape to infinity with probability one. 
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